Optical nonlinearity in Ar and N2 near the ionization threshold 
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We directly measure the nonUnear optical response in argon and nitrogen in a thin gas target 
to laser intensities near the ionization threshold. No instantaneous negative nonlinear refractive 
index is observed, nor is saturation, in contrast with a previous measurement [Loriot et a/.. Opt. 
Express 17, 13429 (2009)] and calculations [Bree et a/., Phys. Rev. Lett. 106, 183902 (2011)]. 
In addition, we are able to cleanly separate the instantaneous and rotational components of the 
nonlinear response in nitrogen. In both Ar and N2, the peak instantaneous index response scales 
linearly with the laser intensity until the point of ionization, whereupon the response turns abruptly 
negative and ^constant, consistent with plasma generation. 



The optical Kerr effect, the intensity-dependent refrac- 
tive index experienced by an optical pulse in a trans- 
parent medium, plays an important role in phenomena 
from nonlinear propagation in optical fibers [1 to mode- 
locking in pulsed lasers [2 to filamentary propagation 
in condensed media and the atmosphere |3j. A recent 
transient birefringence measurement in the components 
of air reported by Loriot et al. [4 purported to show that 
the optical Kerr effect saturates and then becomes neg- 
ative for intensities greater than 26 TW/cm^. A strong 
higher-order Kerr effect, with a crossover from positive 
to negative nonlinear index at intensities well below the 
ionization threshold, would have a huge impact on the 
nonlinear optics of transparent media, and has inspired 
theoretical works predicting plasma-free light filamenta- 
tion and exotic new effects in light propagation [6 . It 
would overturn the picture most have of the mechanism 
behind long-range filamentary propagation of intense ul- 
trashort pulses - as arising from an interplay between 
self-focusing due to the positive optical nonlinearity from 
bound electrons and defocusing due to the plasma gener- 
ated by ionization. The existence of a higher-order Kerr 
effect would also have implications for the general non- 
linear susceptibility in transparent media fP, ^ , including 
harmonic generation [9HTT]. 

Subsequent experimental studies of light filaments [12]- 
fT5] have not supported the higher-order Kerr model, with 
one exception [16 . One measurement [13 found that the 
electron density was two orders of magnitude higher than 
predicted by a calculation including higher-order nonlin- 
earities, but agreed with a simulation based on plasma 
defocusing alone [5 . A physical mechanism for the sat- 
uration and negative response was proposed based on 
the nonlinear response near the threshold of ionization 
pT[ lis] . What is missing from this debate is a direct 
measurement of the nonlinearity that corroborates or re- 
futes the intensity dependence observed by Loriot et al 
Here, we describe such a measurement in Ar and N2 us- 
ing spectral interferometry. We find no saturation and 
no negative instantaneous nonlinear phase, in contrast 
to the original experiment [4 . 

The technique we use, single-shot supercontinuum 



spectral interferometry [19 , provides a single-shot mea- 
surement of the transverse space- and time-dependent 
phase shift of a chirped probe pulse due to a transient 
nonlinearity induced by a short pump pulse. The time 
resolution is given by the probe bandwidth (in this exper- 
iment ^15 fs) and the transverse spatial resolution here is 
3 fim. Loriot et al. used a non-spatially resolved multi- 
shot technique limited in time resolution by the probe 
duration of ~90 fs . They measured the transient bire- 
fringence and inferred the higher-order Kerr coefficients 
from the tensorial symmetry of the nonlinear suscepti- 
bilities x*^^\ X*^^\ [4, 20 . In contrast, we can mea- 
sure parallel and perpendicular components of the non- 
linear response independently. Previously, supercontin- 
uum spectral interferometry was used with 110 fs pump 
pulses to study the nonlinear response of air consituents 
[2T] , and no sign reversal of the nonlinear index was ob- 
served up to intensities where ionization occurs. How- 
ever, the use of a gas cell complicates the interpretation 
of the experiment when the response is highly nonlinear 
[22] . Also, at the high pressures used (> 3 atm) plasma- 
induced refraction limits the peak intensity. Here, we 
measure the nonlinear response of Ar and N2 using a 38 
fs pump pulse and a 2 mm thick gas target. 

Figure [1] shows the experimental setup and data. De- 
tailed descriptions of the experimental technique have 
been given previously [191 [21]. The laser is a 1 kHz repeti- 
tion rate Tiisapphire amplifier producing 38 fs full width 
at half maximum (FWHM), 3.5 mJ pulses centered at 800 
nm. Roughly 700 /iJ of the laser output is used to gen- 
erate supercontinuum covering 640-720 nm in a gas cell 
(not shown) filled with 1-2 atm of Ar; the fundamental is 
rejected using a dichroic mirror. The supercontinuum is 
linearly polarized. Probe and reference pulses, separated 
by 1.4 ps, are generated using a Michelson interferometer 
and chirped so that the group delay dispersion is 1950 fs^, 
and then the beam is spatially filtered with a 100 /im pin- 
hole. The pump power is attenuated using a waveplate 
and thin film polarizers so that the pulse energy is con- 
tinuously adjustable from 5-100 /iJ. The pump beam is 
then expanded with a telescope, and a A/2 waveplate al- 
lows rotation of the pump polarization. The pump and 
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FIG. 1: Experimental apparatus and results, (a) Simplified 
experimental setup diagram showing the pump, probe, and 
reference pulses focused on a flowing gas target in a vacuum 
chamber, (b) Example interferogram. (c) Map of phase versus 
time and transverse dimension for Ar at /peak = 60 TW/cm^, 
extracted as described in [19]. (d) Lineout of the pump spot, 
comparing the signal using the drilled flow tube to the signal 
with a backfilled chamber. The true pump spot is also shown 
for comparison. The dashed lines show the results of propaga- 
tion simulations 22 for the case of a backfilled chamber and 
for an interaction length of 2 mm. (e) The measured peak 
phase shift for Epump || Eprobe as a function of peak intensity 
for Ar (x) and N2 (o). The calculated peak phase shift in Ar 
using the higher-order Kerr coefficients given in [4 is shown 
as a dashed line. 



supercontinuum are combined using a dichroic mirror; at 
this point the pump beam is about 5 times wider than 
the probe/reference beam so that the probe spot overfills 
the pump spot in the interaction region. 

A well-defined interaction length, ideally comparable 
to or shorter than the Rayleigh range of the pump beam, 
minimizes spatial and temporal distortions due to axial 



variation in the pump intensity [22^. To achieve this, we 
use a thin gas target inside a vacuum chamber, shown 
in Fig. [1^. The beams are focused using a lens of focal 
length 40 cm, and enter the vacuum chamber through 
a 5 mm thick fused silica window, propagating 30 cm 
through vacuum to the target. The pump beam waist, 
measured by imaging the focus onto a CCD camera, is 22 
/im FWHM. The target is a copper gas flow tube with a 
flattened section through which a 120 jam diameter hole 
was laser drilled to allow the pump and probe to pass 
through. The space between the inner tube walls is 1 mm, 
with a wall thickness of 0.5 mm. A roots pump keeps the 
background pressure of the chamber at 400 mTorr, and 
the local gas density drops quickly enough away from the 
tube that the probe phase shift is dominated by the gas 
between the holes. The pump beam is rejected at the 
exit of the vacuum chamber using a dichroic mirror that 
is also used as a window. In the data shown here, we 
have subtracted a very small background signal due to 
cross phase modulation in the entrance window. 

The central plane of the gas tube is imaged onto the 
entrance slit of an imaging spectrometer. Interference 
fringes appear in the spectrum because of the time de- 
lay between the probe and reference pulses. An example 
spectral interferogram for Ar at 60 TW/cm^ pump vac- 
uum intensity is shown in Fig. ^p. The intensity values 
given here are calculated from the average power and the 
measured beam profile - we estimate an uncertainty of 
20%. The pump pulse, centered at time t = tpump and 
transverse dimension x = Xpump, causes a phase shift 



in the probe pulse A^(x, t) (too small to produce fringe 
shifts visible by eye in Fig. [i]d). This causes a change 
in the spectral phase and amplitude of the probe beam 
[19]. The spectral phase is found by Fourier analysis of 
the interference fringes, and the change in amplitude is 
also found from the interferogram. The final piece of in- 
formation required is the spectral phase of the reference 
pulse, which is, to an excellent approximation, quadratic 
and proportional to the group delay dispersion [19] . The 
extracted time domain phase shift A0(x, t) of the probe 
is shown in Fig. [l]^. The signal-to-noise ratio is con- 
siderably improved by summing multiple interferograms 
before performing the phase extraction [21] ; in all of the 
data presented here, 300 interferograms were summed at 
each power and polarization. 

The bound electron optical nonlinearity in Ar is in- 
stantaneous to a very good approximation because the 
energy of the lowest electronic excitation is 15 eV, far 
greater than the photon energy 1.5 eV. The ordinary in- 
stantaneous Kerr effect is linear in the intensity /: the 
refractive index is of the form n = no + n2/, where no is 
the index of refraction and n2 is the Kerr coefficient, and 
thus A(j){x^t) (X /pump(^,^)- The measured FWHM in 
Fig. [ij^ of 38 fs matches an autocorrelation measurement 
of the pump pulse. The spatial profile lineout, shown in 
Fig. [iji, agrees well with the pump spot profile. This 
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confirms that the intensity profile in the thin gas target 
is the same as the vacuum profile. The deleterious ef- 
fect of excessive interaction length on the width of the 
response is illustrated in Fig. [T]i. A lineout of the mea- 
sured phase shift along x at t = tpump is shown with the 
chamber backfilled with Ar - note the wider profile com- 
pared to the flow tube case. Also shown are simulations 
of A0(x, tpump) using the beam propagation method [22j 
with the simulation we obtain an effective pump-gas in- 
teraction length of ~2 mm, in agreement with the tube 
geometry. Because of the short interaction length, inten- 
sity clamping pj does not affect the intensity profiles. 

The peak phase shift measured in the experiment is 
plotted as a function of peak intensity in Fig.[l^. We find 
a very linear dependence for both Ar and N2, and nearly 
the same Kerr coefficient, which is consistent with other 
experiments [23|. For peak intensity /peak = 60 TW/cm^, 
we measure a peak phase shift in Ar of 0.059 radians. 
For a medium with an effective interaction length Leff, 
A^peak = 47rLeffn2/peak/A (uotc the extra factor of 2 be- 
cause we measure cross phase modulation). Using the 
literature value for Ar, n2 = 9.8 x 10-^° cm^W [23ll24l 
at 1 atm, and using Leff ~ 2 mm, we estimate an average 
pressure in the interaction region of 0.3 atm. A higher- 
order Kerr effect would add terms of the form n2m/^(^), 
where m > 1 [4 . No negative instantaneous phase is ob- 
served at any intensity in Ar or N2 , nor do we see evidence 
of saturation [TBj , in disagreement with the results of Lo- 
riot et a/. [4]. A simulation of the phase shift expected 
using the coefficients reported in [4] is shown as a dashed 
line in Fig. [l^. The difference is stark and well outside 
any error in our experiment we can conceive of. We have 
also performed the same experiment in Ar using a probe 
pulse whose spectrum overlaps the pump pulse, with or- 
thogonal polarization so that the pump light could be 
rejected by a polarizer before the spectrometer (see the 
appendix). We have studied the possible origin of the re- 
sults obtained in [4 and have found that the interference 
of pump and probe pulses of the same wavelength can 
produce a plasma grating which gives rise to an effective 
birefringence [25 . 

Increasing the pump intensity beyond the level of Fig.[l] 
requires careful consideration of increased supercontin- 
uum generation by the pump itself. Because the pump 
and probe paths are not phase stable with respect to one 
another, spectral fringes between the pump supercontin- 
uum and the reference pulse average out when many in- 
terferograms are summed. So the pump supercontinuum 
does not cause significant data distortion until it satu- 
rates the CCD camera, which occurs at intensities higher 
than 200 TW/cm^, well beyond the ionization threshold 
of -100 TW/cm^ Maps showing A0(x,t) at high 

intensity in Ar are shown in Fig. [2] At high intensi- 
ties we observe an additional response due to the plasma 
generated by ionization [2T. The plasma produces a 
negative index contribution Anpiasma = — ^e/(2A^cr), 



20 
3 40 
I 60 

CO 

o 

^ 80 
100 



0.1 



(a) 




I 



iO.1 



C/i 
CD 

^ 

Q. 

Q)' 

-0.1^ 



200 400 600 
time (fs) 



200 400 600 
time (fs) 



CO 

c 

CO 

TD 
CO 



0^ 



CD 

J5-0.1 



-0.2 



(C) 



15TW/cm^ 

/ I: 1 



120 TW/cm2- 



180 TW/cm2- 



(d) 



15TW/cm^ 



120 TW/cm2 ■ 



180 TW/cm2- 



300 
time (fs) 



600 



300 
time (fs) 



600 



FIG. 2: Experimental data for Ar. (a,c) Epump || Eprobe- 
(b,d) Epump -L Eprobe- (a,b) Extracted phase A(/)(x,t) at 150 
TW/cm^ vacuum pump intensity, showing the positive elec- 
tronic Kerr effect signal at t = tpump and the negative plasma 
signal at later time delays. (c,d) Lineouts A(/)(xpump, t) as 
a function of pump intensity. Within error, the plasma re- 
sponse is constant in time after the pump has passed; we 
observe increased noise at the edges of the time window due 
to the decreased magnitude of the probe/reference spectra on 
the wings. The curves at 120 TW/cm^ and 180 TW/cm^ give 
electron densities Ne — 7x 10^ 
respectively. 



where is the electron density and Ncr is the criti- 
cal density. The plasma densities measured are consis- 
tent with calculations using Ammosov-Delone-Krainov 
(ADK) rates [28 . Plots of A0(xpump,^) as a function 
of pump intensity are shown in Fig. [2j3 (Ep^mp || Eprobe) 
and Fig. [2]i (Epump ^ Eprobe)- Note that the plasma 
contribution seen in Fig. [2] is highly characteristic: (1) 
its onset at higher intensity (180 TW/cm^ compared to 
120 TW/cm^) increasingly dominates the instantaneous 
Kerr response at the back of the pulse; the residual pos- 
itive Kerr peak appears to move forward in time, (2) 
unlike the Kerr response, the plasma-induced phase shift 
is probe polarization independent, and (3) after genera- 
tion, the plasma response is long-lived on the time scale 
of this measurement, owing to recombination timescales 
of order ~100 ps. 

In N2, the optical Kerr response has an additional con- 
tribution from the transient alignment of the molecules 
in the strong optical field [HI [H [30] . Results for N2 are 
shown in Fig. [3] A(/)(x, t) is shown at low pump intensity 
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for parallel and perpendicular polarization in Fig. [3^b. 
The index change is An{t) = n2/(t)+ R{t')I{t-t')dt' , 
where R{t) is a response function that depends on proper- 
ties of the rotational levels and the nuclear spin statistics 
[29l|30]. In N2 the rotational response peaks about 80 fs 
after the pump pulse arrives, as can be seen in Fig. [3^. 
The ratio of the instantaneous Kerr effect for parallel to 
perpendicular polarization is 3:1 in an isotropic medium. 
For the rotational component, the ratio is 2 : —1. The 
different symmetry properties allow the clean separation 
of the two contributions, as shown in Fig. [Sj). Previous 
measurements using this technique [21] |30] were unable 
to resolve the two contributions, but here we can owing 
to the shorter pump pulse. To our knowledge, this is the 
first direct observation of the relative contributions of the 
instantaneous and rotational components of the Kerr ef- 
fect in N2. At higher intensities plasma is observed, as 
shown in plots of A^(xpump,^) in Fig. [sjie. 

In summary, we have performed direct measurements 
of the optical Kerr effect in Ar and N2 using single- 
shot supercontinuum spectral interferometry. In N2, we 
are able to distinguish between instantaneous and rota- 
tional components of the nonlinearity, and the polariza- 
tion dependence is consistent with theory. We observe 
the usual optical Kerr effect, linear in the intensity, as 
well as the onset of plasma, but no higher-order instan- 
taneous nonlinearities effecting either saturation or nega- 
tive response. The fact that the pump-induced response 
appears to be linear in the intensity until the point of ion- 
ization is a reflection of the latter's extremely nonlinear 
onset. At least for 38 fs pulses, there is no practical dis- 
tinction between plasma and special atomic states with 
negative polarizability. Finally, these results disprove the 
idea that higher-order instantaneous nonlinearities are 
important in nonlinear optics in gases at high intensities. 
Our results strongly confirm the long-standing concep- 
tual picture [3 that short pulse filamentation in gases 
arises from the interplay between nonlinear self-focusing 
from bound electron nonlinearities and defocusing due to 
plasma generation. 
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Appendix 

Figure |4] shows experimental data in Ar using probe 
and reference pulses whose power spectra are identical 
to the pump pulse. In this case the probe and reference 
pulses are not supercontinua, but rather chirped versions 
of the pump pulse. To reject most of the pump light 
before the spectrometer, the pump polarization is per- 
pendicular to the probe polarization, and a polarizer is 
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FIG. 3: Experimental data for N2. (a,d) Epun 
(b,e) Epump ^ Eprobe- (a,b) Extracted phase A(/)(x,t) at 60 
TW/cm^ pump intensity, showing the positive instantaneous 
Kerr effect (coincident with the pump) and the rotational 
response at later time delays, (c) Lineouts A(/)(xpump, t) for 
the data shown in (a),(b). The decomposition of the signal 
into instantaneous and rotational components, as described 
in the text, is also shown. (d,e) Lineouts A(/)(xpump, t) as a 
function of pump intensity. 



inserted before the spectrometer. The signal from the 
windows is larger with a probe centered at 800 nm be- 
cause the dichroic mirror that is used as an exit window 
when the probe is supercontinuum centered at 680 nm is 
replaced with a fused silica window. The window contri- 
bution has been subtracted in the data shown in Fig. [4j 
Due to the leakage of pump light through the polarizer 
due to its finite extinction ratio, a small fraction of the 
pump light reaches the CCD camera in the spectrometer. 
We attribute the oscillations in this data to the interfer- 
ence of this residual pump light with the reference pulse; 
as with the supercontinuum generated by the pump, this 
interference decreases when interferograms are averaged 



5 



0.04 



0.02 



■D 

03 



CO 

8? -0.02 

CL 



-0.04 




-0.06 I ' ' ' 1 

200 400 600 

time (fs) 

FIG. 4: Experimental data for Ar with probe and reference 
spectra centered at 800 nm and Epump -L Eprobe- The mea- 
sured phase shift A(/)(xpump, t) is shown, extracted after av- 
eraging 1000 interferograms, at four pump intensities: 25 
TW/cm^ (black), 50 TW/cm^ (red) , 75 TW/cm^ (blue), 
and 100 TW/cm^ (green). 

because the pump and probe arms are not phase sta- 
bilized. In this data, 1000 interferograms were averaged 
before the phase shift was extracted. The time resolution 
with the 800 nm probe is ~ 40 fs because of the narrower 
bandwidth. As with the data taken with a supercontin- 
uum probe, a monotonically increasing, positive instan- 
taneous phase shift is seen up to the ionization thresh- 
old. Note that coherent scattering from a plasma grating 
caused by interference between the pump and probe [25] 
is not present in this data because the pump polarization 
is perpendicular to the probe polarization. 
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